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Multifractal properties of the distribution of topological invariants for a 
model of trajectories randomly entangled with a nonsymmetric lattice of obsta- 
cles are investigated. Using the equivalence of the model to random walks on a 
locally nonsymmetric tree, statistical properties of topological invariants, such 
as drift and return probabilities, have been studied by means of a renormaliza- 
tion group (RG) technique. The comparison of the analytical RG-results with 
numerical simulations as well as with the rigorous results of P.Gerl and W.Woess 



[15 1 demonstrates clearly the validity of our approach. It is shown explicitly by 
direct counting for the discrete version of the model and by conformal methods 
for the continuous version that multifractality occurs when local uniformity of 
the phase space (which has an exponentially large number of states) has been 
broken. 



I. INTRODUCTION 

The phenomenon of multifractality consists in a scale dependence of critical exponents. It 
has been widely discussed in the literature for a wide range of issues, such as statistics of 
strange sets IQ, diffusion limited aggregation 0, wavelet transforms 0, conformal invariance 
Q or statistical properties of critical wave functions of massless Dirac fermions in a random 
magnetic field [^|-|^. 

The aim of our work is not only to describe a new model possessing multiscaling dependence, 
but also to show that the phenomenon of multifractality is related to local nonuniformity of the 
exponentially growing ("hyperbolic") underlying "target" phase space, through an example of 
entangled random walks distribution in homotopy classes. Indeed, to our knowledge, almost 
all examples of multifractal behavior for physical or more abstract [|T],^ systems share 

one common feature — all target phase spaces have a noncommutative structure and are locally 
nonuniform. 

We believe that multiscaling is a much more generic physical phenomenon compared with 
uniform scaling, appearing when the phase space of a system possesses a hyperbolic structure 
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with local symmetry breaking. Such perturbation of local symmetry could be either regular 
or random — from our point of view the details of the origin of local nonuniformity play a less 
significant role. 

We discuss below the basic features of multifractality in a locally nonuniform regular hy- 
perbolic phase space. We show in particular that a multifractal behavior is encountered in 
statistical topology in the case of entangled (or knotted) random walks distribution in topolog- 
ical classes. 

The paper is organized as follows. In Section |I| we consider a 2D A^-step random walk in a 
nonsymmetric array of topological obstacles and investigate the multiscaling properties of the 
"target" phase space for a set of specific topological invariants — the "primitive paths". The 
renormalization group computations of mean length of the primitive path, as well as return 
probabilities to the unentangled topological state are developed in Section |ITT[ Section is 
devoted to the application of conformal methods to a geometrical analysis of multifractality in 
locally nonuniform hyperbolic spaces. 



II. MULTIFRACTALITY OF TOPOLOGICAL INVARIANTS FOR RANDOM 
ENTANGLEMENTS IN A LATTICE OF OBSTACLES 

The concept of multifractality has been formulated and clearly explained in the paper . We 
begin by recalling the basic definitions of Renyi spectrum, which will be used in the following. 

Let i^(Cj) be an abstract invariant distribution characterizing the probability of a dynamical 
system to stay in a basin of attraction of some stable configuration Ci {i = 1,2, .. . ,Af). Taking 
a uniform grid parameterized by "balls" of size I, we define the family of fractal dimensions Dq-. 



Da = lim (1) 

^ g - 1 In/ ^ ^ 

As q is varied, different subsets of u'^ associated with different values of q become dominant. 
Let us define the scaling exponent a as follows 

where a can take different values corresponding to different regions of the measure which become 
dominant in Eq.(|^). In particular, it is natural to suggest that i^'^{Ci) can be rewritten as 
follows: 



AT r. 
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where p{a) is the probabihty to have the value a lying in a small "window" [a', a' + Aa'] and 
/(a) is a continuous function which has sense of fractal dimension of the subset characterized 
by the value a. 

Supposing p{a) > one can approximately evaluate the last expression via the saddle-point 
method. Thus, one gets (see, for example, [||]): 

-J{a) < 



what together with (|I]) leads to the following equations 

r(g) = qa{q) - f[a{q)\ 

where r(g) = (g — 1)-Dg. Hence, the exponents r(g) and /[a(g')] are related via Legendre 
transform. For further details and more advanced mathematical analysis, the reader is refered 

to m. 



A. 2D topological systems and their relation to hyperbolic geometry 



Topological constraints essentially modify physical properties of the broad class of statis- 
tical systems composed of chain-like objects. It should be stressed that topological problems 
are widely investigated in connection with quantum field and string theories, 2D gravitation, 
statistics of vortices in superconductors, quantum Hall effect, thermodynamic properties of 
entangled polymers etc. Modern methods of theoretical physics allow us to describe rather 
comprehensively the effects of nonabelian statistics on the physical behavior of some systems. 
However the following question remains still obscure: what are the fractal (and as it is shown 
below, multifractal) properties of the distribution function of topological invariants, character- 
izing the homotopy states of a statistical system with topological constraints? We investigate 
this problem in the framework of the model "Random Walk in an Array of Obstacles" (RWAO). 

The RWAO-model can be regarded as physically clear and as a very representative image for 
systems of fluctuating chain-like objects with a full range of nonabelian topological properties. 
This model is formulated as follows: suppose that a random walk of N steps of length a takes 
place on a plane between obstacles which form a simple 2D rectangular lattice with unit cell of 
We assume that the random walk cannot cross ("pass through") any obstacles. 
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(a) (b) 

FIG. 1. Random walk in the two-dimensional rectangular lattice of obstacles. 



It is convenient to begin with the lattice realization of the RWAO-model. In this case the 
random path can be represented as a A^-step random walk in a square lattice of size a x a 
{a < Cy < Cx) — see fig.|]. 




(a) (b) 

FIG. 2. Random path for a = Cx = Cy-. (a) in the 2D lattice of obstacles; (b) in the covering space 
(on the Cayley tree). 



Cy a lattice 



It had been shown previously (see, for example |0,|TT1) that for a = Cx 
random walk in the presence of a regular array of obstacles (punctures) on the dual lattice 

is topologicaly equivalent to a free random walk on a graph — a Cayley tree with branching 
number 2; = 4 (see Fig.^. An outline of the derivation of this result is as follows. The different 
topological states of our problem coincide with the elements of the homotopy group of the 
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multi-punctured plane, which is the free group Too generated by a countable set of elements. 
The translational invariance allows to consider a local basis and therefore to study the factored 
group Foo/Z^ = Tz/2, where Tz/2 is a free group with z/2 generators whose Cayley graph is 
precisely a 2;-branching tree. 

The relation between Cayley trees and hyperbolic geometry is discussed in details in Sec- 
tion |T|. Intuitively such a relation could be understood as follows. The Cayley tree can be 
isometrically embedded in the hyperbolic plane Ti. (surface of constant negative curvature). 
The group Tz/2 is one of the discrete subgroups of the group of motion of the hyperbolic plane 
TC = SL{2,R)/ S0{2), therefore the Cayley tree can be considered as a particular discrete 
realization of the hyperbolic plane. 

Returning to the RWAO-model, we conclude that each trajectory in the lattice of obstacles 
can be lifted to a path in the "universal covering space" i.e to a path on the z-branching 
Cayley tree. The geodesic on the Cayley graph, i.e the shortest trajectory along the graph 
which connects ends of the path, plays the role of a complete topological invariant for the 
original trajectory in the lattice of obstacles. For example, the random walk in the lattice of 
obstacles is closed and contractible to a point (i.e. is not entangled with the array of obstacles) 
if and only if the geodesic length between the ends of the trajectory on the Cayley graph is 
zero. Hence, this geodesic length can be regarded as a topological invariant, which preserves 
the main nonabelian features of the considered problem. 

We would like to stress two facts concerning our model: (i) The exact configuration of a 
geodesic is a complete topological invariant, while its length k is only a partial topological 
invariant (except the case k = 0); (ii) Geodesies have a clear geometrical interpretation, having 
sense of a bar (or "primitive") path which remains after deleting all even times folded parts 
of a random trajectory in the lattice of obstacles. The concept of "primitive path" has been 
repeatedly used in statistical physics of polymers, leading to a successful classification of the 
topological states of chain-like molecules in various topological problems |]TU|-|T^. 

Even if many aspects of statistics of random walks in fixed lattices of obstacles have been 
well understood (see, for example and references therein), the set of problems dealing with 
the investigation of fractal properties of the distribution of topological invariants in the RWAO- 
model are practically out of discussion. Thus we devote the next Section to the study of fractal 
and multifractal structures of the measure on the set of primitive paths in the RWAO-model 
for a Cy < Cx- 



B. Multifractality of the measure on the set of primitive paths on a nonsymmetric 

Cayley tree 

The classification of different topological states of a A^-step random walk in a rectangular 
lattice of obstacles in the case a -C c„ < c^; turns out to be a more difficult and more rich 
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problem than in the case a = Cy = Cx discussed above. However, after a proper rescahng, the 
mapping of a random walk in the rectangular array of obstacles to a random walk on a Cayley 
tree can be explored again. To proceed we should solve two auxiliary problems. First of all we 
consider a random walk inside the elementary rectangular cell of the lattice of obstacles. Let 
us compute: 

(i) The "waiting time", i.e the average number of steps (t) which a t-step random walk 
spends within the rectangle of size x Cy] 

(ii) The ratio of the "escape probabilities" and Py through the corresponding sides and 
Cy for a random walk staying till time t within the elementary cell. 

The desired quantities can be easily computed from the distribution function P(xo, yo, x, y, t) 
which gives the probability to find the t-step random walk with initial (xo,?/o) and final {x,y) 
points within the rectangle of size x Cy. The function P{x,y,t) in the continuous approxi- 
mation (a ^ 0; t — » oo; at = const) is the solution of the following boundary problem 



P(0, y, t) = P{cx, y, t) = P(x, 0, t) = P(x, Cy, t) = 
. Pix,y,0) = 6{xo,yo) 



(3) 



where a is the length of the effective step of the random walk and the value ^ has sense of a 
diffusion constant. 

The solution of Eqs.(|]) reads 



P{xo,yo,x,y,t) 



A 00 00 

m^=l my=l 



2 2 

T a 
4 



'X '-y 



Tim^XQ TirriyyQ Tim^x 7rm„?/ 
sin sin — sin sin — 



(4) 



The "waiting time" (t) can be written now as follows 

1 rc: 



it) 



Cx Cy JO 



rcy pCx rcy roo 

dxo dyo dx dy dt P{xo,yo,x,y,t) 
Jo Jo Jo Jo 



(5) 



while the ratio Px/Py can be computed straightforwardly via the relation: 



r^^x f^y r^x 

dxo dyo dx P{xo,yo,x,y,t) 
Px Jo Jo Jo 



y={a,Cy-a} 



Py 



Cx rCy l-Cy 

dxo / dyo / dy P{xo,yo,x,y,t) 
Jo Jo 



x={a,Cx — a} 



(6) 
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In the "ground state dominance" approximation we truncate the sum (^) at rrix = rriy = 1 
and get the following approximate expressions: 

4*^ n 



TT^ {cl + c2) ' Py C: 



In the symmetric case {c^ = Cy = c) Eq.(^ gives (t) = |^ and Px/Py = 1, as it should be for 
the square lattice of obstacles. 

Now the distribution function of the primitive paths for the RWAO model can be obtained 
via lifting this topological problem to the problem of directed random walksQ on the 4-branching 
Cayley tree, where the random walk on the Cayley tree is defined as follows: 

(a) The total number of steps N on the Cayley tree is 



N 



(t) 44 cl a 



2 



-x ^y 

(the value (t) has been computed in (|^)). 

(b) The distance (or "level" k) on the Cayley tree is defined as the number of steps of the 
shortest path between two points on the tree. Each vertex of the Cayley tree has 4 branches; the 
steps along two of them carry a Boltzmann weight 1, while the steps along the two remaining 
ones carry a Boltzmann weight j3 as it is shown in Fig.|. The value of /? is fixed by Eq.(0), 
which yields 

/3 = ^ = 4 (8) 
Py S 

The ultrametric structure of the topological phase space, i.e. of the Cayley tree 7(/3), allows 
us to use the results of paper ^ for investigating multicritical properties of the measure of 
all primitive (directed) paths of k steps along the graph 7(/?) with nonsymmetric weights 1 
and P (see Fig.^. A rigorous mathematical description of such weighted paths on trees (called 



cascades) can be found in ||25|, where the authors derive multifractal spectra, but for different 
distributions of weights. 



Recall that by definition the primitive path is the geodesic distance and therefore cannot have two 
successive opposite steps. 
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1/31 /31/3 1/31 /31/3 1/31 f3 1 (3 1/3 1 /3 1 /3 1/3 1 /3 1 /3 1/3 1 /5 1 /3 

FIG. 3. 4-branching Cayley tree with different transition probabilities along branches. 

We construct the partition function k) which counts properly the weighted number of 
all 4 X 3*^"^ different fc-step primitive paths on the graph 7(/5). 

Define two partition functions a/^ and bk of fc-step paths, whose last steps carry the weights 
1 and (3 correspondingly. These functions satisfy the recursion relations for A; > 1: 



bk+i = 2pak + pbk 



{k>l] 



(9) 



with the following initial conditions at A; = 1: 



ai = 2 
bi = 2(3 



(10) 



Combining (P) and ([l0|) we arrive at the following 2-step recursion relation for the function 



ak+2 = (1 + /5) ttfc+i + 3/3 ttk {k > 1) 
ai = 2 {k = l) 

a2 = 2 + 4/3 {k = 2) 



whose solution is 



ak 



0,2 " 01A2 , o-iAi — a2 



Ai — A2 



Ar^ + 



Ai — A2 



(12) 



where 



Ai,2 = -(l + /5±V(l+/5)' + 12/5 



(13) 
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Taking into account that bk is given by the same recursion relation as but with the initial 
values bi = 2(3 and 62 = 2/5^ + 4/5, we get the following expression for the partition function 
Vt{(3, k) =ak + bk. 



2(1 +Af3 + (3^) - 2(1 + /?)A2 ^k-i , 2(1 + /3)Ai - 2(1 + A(3 + (3^) 



Ai — A2 



Ai — A2 



(14) 



The partition function VL{(3, k) contains all necessary information about the multifractal 
behavior. Following Eqs.(|I])-(|D, we associate the set of stable configurations {Cj} with the set 
of Af{k) = 4 X 3'^"^ vertices of level k. Hence, we define 



(15) 



Boltzmann weight = 0.001 



Boltzmann weight = 0.01 



0.00375 - 



0.0037 - 



0.00365 



0.03482 



0.0348 



0.03478 



0.03476 



0.03474 



0.03472 



0.0036 



I , ; , ; , ; , ; , 1 0.0347 
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0.25703 



0.25702 



0.25701 



0.257 - 



0.25699 - 



Boltzmann weight = 0.1 



0.711947 



0.711942 



0.711937 



Boltzmann weight = 0.5 



0.25698 ' ' ' ' ' ' ' ' ' ' ' 0.711932 

-25 -15 -5 5 15 25 -25 -15 -5 5 15 25 



FIG. 4. Dependence Dq[q) for different values of Boltzmann weight (3. 



Taking into account that the uniform grid has resolution l{k) = l/7V(fc) for > 1 and using 
Eq.(|I]), we obtain 



T{q) = hm 



\nl{k) 



(16) 



which allows to determine the generalized Hausdorff dimension Dg via the relation 

Dg = Tiq)/iq-l) 



(17) 



The corresponding plots of the functions Dq{q) for different values oi (3 = 
{0.001; 0.01; 0.1; 0.5} are shown in Fig.|| (the numerical computations of Eqs. ([T6|)-(p!7D are 
carried out for k = 100 000). The fact that Dq{q) depends on q clearly demonstrate the multi- 
fractal behavior. 



III. RANDOM WALK ON A NONSYMMETRIC CAYLEY TREE 



A. Master equation 



Consider a random walk on a 4-branching Cayley tree and investigate the distribution 
P{k, N) giving the probability for a A^-step random walk starting at the origin of the tree to 
have a primitive (shortest) path between ends of length k. The random walk is defined as 
follows: at each vertex of the Cayley tree the probability of a step along two of the branches is 
Px, and is Py along the two others; p^ and Py satisfy the conservation condition + 2py = 1. 
Using Eq.(|D, the following expressions hold: 



Px 



Py 



/3 



2(1 + /?) 
1 

WW) 



(18) 



The symmetric case (3=1 (which gives p^ = Py = 1/4) has already been studied and an 
exact expression for P{k,N) has been derived in [|l^. The rigorous mathematical description 
of random walks on graphs can be found in . Importance of spherical symmetry (i-e the fact 
that all vertices of a given level are strictly equivalent) is discussed in Another example 
of nonsymmetric model on a tree (case of randomly distributed transition probabilities, the so- 
called RWRE model) is described in To our knowledge, the solution for the nonsymmetric 



random walk which we defined above is known only for k fixed and ^ 1 |jT5[ . Here we consider 
the case k ^ 1, N ^ 1, and in particular we study the distribution in the neighborhood of 
the maximum. Breaking the symmetry by taking (3^1 affects strongly the structure of the 
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problem, since then the phase space becomes locally nonuniform: we have now vertices of 
two different kinds, x and j/, depending on whether the step toward the root of the Cayley 
tree occurs with probability or py. In order to obtain a master equation for P{k,N), we 
introduce the new variables Lx{k, N) and Ly{k, N), which define the probabilities to be at the 
level in a vertex x or y after steps. We recursively define the same way the probabilities 
Lai...a„{k, N), (oj = {x, y}) to be at level k in a. vertex such that the sequence of vertices toward 
the root of the tree is ai...a„. One can see that the recursion depends on the total "history" 
till the root point, what makes the problem nonlocal. The master equation for the distribution 
function P{k, N) 

P{k,N + l) = {2p,+py)Ly{k-l,N) + {2py+p,)L,{k-l,N) + 

PyLy{k+l,N)+P^L^{k+l,N) 

is coupled to the hierarchical set of functions {L^, Ly] L^^, L^y, Ly^, Lyy] La^^ a^} which satisfy 
the following recursion relation 

La,...aAk, N + 1) = (2 - 5a„a,)Pa,La,...a^ik - 1, N)+ ^ ^^Q) 

xai...a„ 

{k + l,N) +pyL 

ya\...an 

where ai...a„ cover all sequences of any lengths (< k) in x and y. In order to close this infinite 
system at an arbitrary order no we make the following assumption: for any n < Hq we have 



P{k, N) 



' aa^...a„ (21) 

k^riQ 



with aa-^...a„ constant. 

Using the approximation (|2T|) we rewrite ([T9|)-(^0D for large k and in terms of the function 
P{k,N) and constants aai...a„ (0 < n < no). Taking into account that 

Lia\...anX ~l~ n^y -^ai...an 

we arrive at 2"""^ independent recursion relations for one and the same function P{k, N), with 
2"o _ I independent unknown constants «ai...a„Q- In order to make this system self-consistent, 
one has to identify coefficients entering in different equations, what yields 2"^° — 2 compatibility 
relations for the constants aai...a„Q, and the system is still open. This fact means that all scales 
are involved and the evolution of La^ a^ depends on i^ai...a„+i, the evolution of La^ a^+i depends 
on Laj...a„^2 and so on. At each scale we need informations about larger scales. This kind of 
scaling problem naturally suggests to use a renormalization group approach, which is developed 
in the next Section. 

To begin with the renormalization procedure, we need to estimate the values of the constants 
Q^ai...a„o for the first (i.e. the smallest) scale. Let us denote 

= a 
ay = 1 — a 
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and define a^x-, o^xy, ctyy, ctyx as follows: 



a 



xy 



yx 



(1 - Vx)a 
{l-vy){l-a) 



Now we set 



Px(^xa\...an Py(^yai...a„ (^Px(^ ~^ Py(,^ 



a] ]a 



ai...a„ 



(22) 



what means that we neglect the correlations between the constants aai...a„ and aa2...a„ at dif- 
ferent scales. As it is shown in the next Section, the renormalization group approach allows us 
to get rid of the approximation (p2D. 

With (|22| ) one can obtain the following generic master equation 

Pai CX.a2...an 



(^ai...a„ 



(2 - 6a,,a,)P{k -l,N) + {apx + (1 - a)py)P{k + 1, N) (23) 



where ai...a„ again cover all possible sequences in x and y. We have now 2"° — 1 unknown 
quantities with 2"" — 1 compatibility relations (P5|), what makes the system ( ^5] ) closed. 

For illustration, we derive the solution for uq = 2: 
P 



(24) 



P{k, N+l) = ^P{k - 1, iV) + [apx + (1 - a)py]P{k + 1, N) 



P{k,N+l] 
P{k,N+l] 
P{k,N+l] 



2p^{l - a] 
Py 



P{k - 1, iV) + (^apx + (1 - a)py)P{k + 1, iV) 



^P(A; - 1, iV) + (ap, + (1 - a)py)P{k + 1, iV) 



2pya 



{l-a){l-vy) 



Note that ( P^ ) displays clearly a Z2 symmetry: Px ^ Py, ct 
conditions for system (p^ read: 

2pa;(l - a) 



P{k - 1, iV) + (ap, + (1 - a)py)P{k + 1, iV) 

—a, Vx Vy. Compatibility 

2pya 



Vx Vy a{l-Vx) {l-a){l-Vy) 



(25) 



which finally gives 



a 



Vx 



-1 - + ^/TTTWTW 



2(1-/?) 



a 



2- a 
1 — a 
1 + a 



(26) 
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(27) 



As it has been said above, without ( p2D the system (p^)-pO|) is open, giving a single equation 
for the unknown function P{k, N) depending on the unknown parameter a: 

P{k,N + l) = (^{2p,+py)(l-a) + {2py+p.,)a)P{k-l,N) + 
{py{l-a) +p,a)P{k + l,N) 

Eg . (pTl) describes a ID diffusion process with a drift 

^ = k = 2apy + 2{l-a)p^ (2^ 



and a dispersion 

5 = = 1 - ^{apy + (1 - aK)' (29) 



which provides for A; ^ 1 and ^ 1 the usual Gaussian distribution with nonzero mean (see 
||1U||). The value of a obtained in ( PBj ) using the approximation (^) gives a fair estimate of the 
drift compared with the numerical simulations, as it is shown in Fig.^. 



B. Real space renormalization 



In order to improve the results obtained above, we recover the information lost in the 
approximation (^2j) and take into account "interactions" between different scales. Namely, we 
follow the renormalization flow of the parameter a{l) at a scale I supposing that a new effective 
step is a composition of 2' initial lattice steps. Let us define: 



• the probability /a(0 of going forth (with respect to the location of the root point of the 
Cayley tree) from a vertex of kind a; 

• the probability ba{l) of going back (towards the root point of the Cayley tree) from a 
vertex of kind a; 

• the probability a{l) of being at a vertex of kind x; 

• the conditional probability Wa{l) to reach a vertex of kind a starting from a vertex of kind 
a under the condition that the step is forth; 

• the conditional probability Va{l) to reach a vertex of kind a starting from a vertex of kind 
a under the condition that the step is back; 

• the effective length d{l) of a composite step. 
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Then the drift k{l) at scale / is given by (compare with |28|) : 

k{l) = d{l) \a{l)(m - bM + (l - a{l))(fy{l) - by{l) 



(30) 



We say that the problem is scale-independent if the flow k{l) is invariant under the decima- 
tion procedure, i.e. with respect to the renormalization group. We compute the flow counting 
the appropriate combinations of two steps, depending on the variable considered: 

Wa{l + 1) -- 
Vail + I) -- 



fa{i + r 

bail + I] 

dil + 1) 
a(/ + 1) 



1-Wail)){l-W,il))+Wlil) 
1-Vail)){l-V,il))+Vlil) 
fail) [Wail) fail) + Wail)) Ml) 



Cail) 

bail)[vail)bail) + {l-Vail))b-ail) 



(31) 



Cail) 

dil) [a(/)c,(/) + (l-a(/))c,(/) 
kil) [ail)w^il) + (l - ail)) (l - Wyil) 



+ 



1 



kil)) [a{l)v^{l) + {l-ail))(l-Vyil) 



where d = x when a = y (and a = y when a = x) and the value Ca(/) ensures the conservation 
condition /a(/ + 1) +6a(/ + 1) = 1 because we do not consider the combinations of two successive 
steps in opposite directions. 

The transformation of a in (|3lD needs some explanations. We consider the drift fc(/) as 
a probability to make a (composite) step forward. The equation for a is given by counting 
the different ways of getting to a vertex of kind x. One can check that /c(/) given by (|30D 
remains invariant under such transformation, what is considered as a verification of the scale 
independence (i.e. of renormalizability) . 

Following the standard procedure, we find the fixed points for the flow of ail). First of all we 
realize that the recursion equations for Wail) and Vail) can be solved independently, providing 
a continuous set of fixed points: w° = 1 — w° and = 1 — Vy. Using the initial conditions (^61) 
for Vail) and deriving straightforwardly the absent initial conditions for Wail), we get 



Vxi^) = 

Vyil) = Vy 

w^il) = 

Wyil) = Wy 



Px 



Px + 
Py 

Py + 2pa; 



(32) 
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(we recall that these values are obtained by taking into account the elementary correlations for 
two successive steps). 

With the initial conditions (|32|) we find the following renormalized values f ° and at the 
fixed point 

1 



= = lim v^{l) = 1 - lim Vy{l) 

J— >oo +00 2 

w° = w\p) = lim w^{l) = 1 - lim Wy{l) = ^ 

(— >oo l—*oo 2 



n=l 



n=l 



where /'•"^(x) is the n^^ iteration of the function 

f{x) =x^ -2x + 2 

We then obtain successively all renormalized values at the fixed point 



J a 


= 1 


IP 


= 




= F 



a 



where the invariance of the drift k is taken into account: 



(33) 



(34) 
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FIG. 5. The values of a and k compared to renormalized quantities and numerical simulations. 
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In figure |^ we compare the theoretical results with numerical simulations. It is worth 
mentioning the efficiency of the renormalization group method, which yields a solution in very 
good agreement with numerical simulations in a broad interval of values (3. 

In addition we compare our results with the exact expression obtained by P.Gerl and 



W.Woess in [|T5[ for the probability P(0, A^) to return to the origin after random steps 
on the nonsymmetric Cayley tree. This distribution function P(0, A^) reads 



P(0, N) oc /i^ 



(35) 



with 



min|-t + y't2 + 4p2 + y't2 + |^ > gj 



(36) 



Let us assert now without justification that Eq.(|27|) (which is actually written for k ^ 1 
and ^ 1) is valid for any values of k and A^. The initial conditions for the recursion relation 
(B7|) are as follows 



(37) 



P(0, N+l) = {2ap, + 2(1 - a)py)P{l, N) 
P{k,0) = 6k,o 



One can notice that Eq.(p7|) completed with the conditions (pT]) can be viewed as a master 
equation for a symmetric random walk on a Cayley tree with effective branching z continuously 
depending on /3: 



ap^ + {1- a)py 



(38) 



Hence, we conclude that our problem becomes equivalent to a symmetric random walk on a 
2; (/?) -branching tree. For A; = the solution, given in [|1^ 



IS 



P(0,iV) oc 



2V^(/3) - 1 



N 



(39) 



This provides the same form as the exact solution (pSj). It has been checked numerically that 
for P G IR^ the discrepancy between (BH) and (B^) is as follows 



< 0.02 



Thus, we believe that our self-consistent RG-approach to statistics of random walks on non- 
symmetric trees can be extended with sufficient accuracy to all values of k. 
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IV. MULTIFRACTALITY AND LOCALLY NONUNIFORM CURVATURE OF 

RIEMANN SURFACES 



We have claimed in Sections |I|-|I| that local nonuniformity and the exponentially growing 
structure of the phase space of statistical systems generates a multiscaling behavior of the cor- 
responding partition functions. The aim of the present Section is to bring geometric arguments 
to support our claim by introducing a different approach of the RWAO model. The differences 
between the approach considered in this Section and the one discussed in Section ^ are as 
follows: 

• We consider a continuous model of random walk topologically entangled with either a 
symmetric or a nonsymmetric triangular lattice of obstacles on the plane. 

• We pursue the goal to construct explicitly the metric structure of the topological phase 
space via conformal methods and to relate directly the nonuniform fractal relief of the 
topological phase space to the multifractal properties of the distribution function of topo- 
logical invariants for the given model. 

Consider a random walk in a regular array of topological obstacles on the plane. As in the 
discrete case we can split the distribution function of all A^-step paths with fixed positions of 
end points into different topological (homotopy) classes. We characterize each topological class 
by a topological invariant similar to the "primitive path" defined in Section |I[ Introducing 
complex coordinates z = x + iy on the plane, we use conformal methods which provide an 
efficient tool for investigating multifractal properties of the distribution function of random 
trajectories in homotopy classes. 

Let us stress that explicit expressions are constructed so far for triangular lattices of obstacles 
only. That is why we replace the investigation of the rectangular lattices discussed in Sections 
by the consideration of the triangular ones. Moreover, for triangular lattices a continuous 
symmetry parameter (such as (3 = c^/Cy in case of rectangular lattices) does not exist and only 
the triangles with angles (vr/S, 7r/3, 7r/3), (7r/2, 7r/4, 7r/4), (7r/2, 7r/6, 7r/3) are available — only 
such triangles tessellate the whole plane z. In spite of the mentioned restrictions, the study of 
these cases enables us to figure out an origin of multifractality coming from the metric structure 
of the topological phase space. 

Suppose that the topological obstacles form a periodic lattice in the z-plane. Let the funda- 
mental domain of this lattice be the triangle ABC with angles either (vr/S, tt/S, vr/S) (symmetric 
case) or (7r/2, 7r/6, 7r/3) (nonsymmetric case). The conformal mapping z{() establishes a one- 
to-one correspondence between a given fundamental domain ABC of the lattice of obstacles 
in the z-plane with a zero-angled triangle ABC lying in the upper half-plane r] > of the 
plane C = + ^-^id having corners on the real axis rj = 0. To avoid possible misunderstand- 
ings let us point out that such transform is conformal everywhere except at corner (branching) 
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points — see, for example [|T^ . Consider now the tessellation of the 2;-plane by means of consec- 
utive reflections of the domain ABC with respect to its sides, and the corresponding reflections 
(inversions) of the domain ABC in the C^plane. Few first generations are shown in Fig.^ The 
obtained upper half-plane Im^ > has a "lacunary" structure and represents the topological 
phase space of the trajectories entangled with the lattice of obstacles. The details of such a 
construction as well as a discussion of the topological features of the conformal mapping z{() 
in the symmetric case can be found in |]13|. We recall the basic properties of the transform z{() 
related to our investigation of multifractality. 




FIG. 6. Conformal mapping of the complex plane z to the "lacunary" upper half-plane Im(^ > 
endowed with a Poincare metric. The zero-angled triangle ABC on corresponds to the triangle ABC 
with the angles (7r/2, 7r/6, 7r/3) on z. 



The topological state of a trajectory C in the lattice of obstacles can be characterized as 
follows. 



• Perform the conformal mappings Zs{C) (or -2ns(C)) of the plane z with symmetric (or 
nonsymmetric) triangular lattice of obstacles to the upper half-plane Im^ > 0, playing 
the role of the topological phase space of the given model. 

• Connect by nodes the centers of neighboring curvilinear triangles in the upper half-plane 
lm( > and raise a graph 7s (or 7ns) (which is, as shown below an isometric Cayley tree 
embedded in the Poincare plane). 

• Find the image of the path C in the "covering space" lm( > and define the shortest 
(primitive) path connecting the centers of the curvilinear triangles where the ends of the 
path C are located. The configuration of this primitive path projected to the Cayley tree 
7s (or 7ns) plays the role of topological invariant for the model under consideration. 

The Cayley trees 7s,ns have the same topological content as the one described in the Section 
H, but here we determine the Boltzmann weights ^1,^2,^3 associated with passages between 
neighboring vertices (see Fig.0) directly from the metric properties of the topological phase 
space obtained via the conformal mappings ^s,ns(C)- 
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FIG. 7. Nonsymmetric 3-branching Cayley tree. 



It is well known that random walks are conformally invariant; in other words the diffusion 
equation on the plane z preserves its structure under a conformal transform, but the diffusion 
coefficient can become space-dependent |18|. Namely, under the conformal transform the 



Laplace operator = is transformed in the following way 



' (40) 



dzd^ |Z'(C)|2 d(dC 

Before discussing the properties of the Jacobians for the symmetric and nonsym- 

metric transforms, it is more convenient to set up the following geometrical context. The 
connection between Cayley trees and surfaces of constant negative curvature has already been 
pointed out [0, mostly through volume growth considerations. Therefore it becomes more 



natural to regard the upper half-plane lm( > as the standard realization of the hyperbolic 
2-space (surface of constant negative curvature R, with here arbitrarily R = —2), that is to 
consider the following metric: 

= ^{de + dv') (41) 

Let us rewrite the Laplace operator 

'^';z = m,v)v' (3^ + 33) (42) 




dzdz 

where the value -D(^, rj) = -D(C) can be interpreted as the normalized space-dependent diffusion 
coefficient on the Poincare upper half-plane: 

^(0 = 1 Lm2 (43) 
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The methods providing the conformal transform ^;s(C) fo^' the symmetric triangle with angles 
(tt/S, 7r/3, 7r/3) have been discussed in details in |]T^. The generalization of these results to the 
conformal transform 2;ns(C) fo^' the nonsymmetric triangle with angles (7r/2, vr/G, vr/S) is very 
straightforward. We here expose the Jacobians of those conformal mappings without derivation: 



l<(C)f 



I4(c)r 



52 I 1 
12' 3 



^0(0, e 



8/3 



8/3 



^2(0, e^'^^)' ^3(0,6^'^^)'^' 



(44) 



where 9[{x\---) = {i = 0, ...,3) are the standard definitions of Jacobi 

elliptic functions |19[. 

Combining ( ^3] ) and (Q) we define the effective inverse diffusion coefficients in symmetric 
(DjT^) and in nonsymmetric {D:^J^) cases: 



(45) 



The corresponding 3D plots of the reliefs ^(^, 77) and ^{$,, rj) are shown in Fig.[ 



symmetric case 



nonsymmetric case 



-log y 2 




-log y 2 




FIG. 8. Profiles of the surfaces ^{C^v) ^iid D^^{^,r]) where ^ = x, rj = — logy. First generation 
of horo cycles are shown by the lines. 



The functions D^^{() and D~^{() are considered as quantitative indicators of the topological 
structure of the phase spaces; in particular a Cayley tree can be isometrically embedded in the 
surface D~^{(). It can be shown that the images of the centers of the triangles of the symmetric 
lattice in the z-plane correspond to the local maxima of the surface D~^{() in the C^plane. We 
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define the vertices of tlie embedded tree as tliose maxima. Tlie links connecting neigfiboring 
vertices are defined in the next paragraph. 

Let us define the horocycles which correspond to repeating sequences of weights in Fig.|^ with 
minimal periods. There are only three such sequences: /5i/52/3i/32---, PiPsPiPa--- (^2[^^(^2[^z---- 
The horocycles are images (analytically known) of certain circles of the 2;-plane. They proved 
to be a convenient tool for a constructive description of the trajectories in the z-plane starting 
from the trajectories in the covering space C,. 

The first generation of horocycles (closest to the root point of the Cayley tree) is shown in 
Fig.^. Let us consider the symmetric case. Following a given horocycle we follow a ridge of the 
surface, and we pass through certain maxima of this surface (that is through certain vertices 
of the tree). We therefore define locally the links of the tree as the set of ridges connecting 
neighboring maxima of D~^{^, if). We recall that the ridge of the surface can be defined as the 
set of points where the gradient of the function D~^{^,7]) is minimal along its isoline. Even if 
this gives a proper definition of the tree, extracting a direct parametrization is difficult, that is 
why henceforth we will approximate the tree by arcs of horocycles. 

To give a quantitative formulation of the local definition of the embedded Cayley tree, we 
consider the path integral formulation of the problem on the (-plane. Define the Lagrangian 
C oc D~^(C)C^ of a free particle moving with the diffusion coefficient -Ds(C) the space (. 
Following the canonical procedure and minimizing the corresponding action pO[, we get the 
equations of motion in the effective potential U = ln[rj'^Ds): 

Qi = {qjdjU)qi - ^qjqjdiU (46) 

where qi = and q2 = f]. Even if Eg. (^6]) is nonlinear with a friction term, one can show that 
the trajectory of extremal action between the centers of two neighboring triangles follow the 
ridge of the surface -Ds(C)- 

It is noteworthy that obtaining an analytical support of Cayley graphs is of great importance, 
since those graphs clearly display ultrametric properties and have connections to p-adic surfaces 
2l| . The detailed study of metric properties of the functions D~^{Q and D~^{() is left for 



separate publication. 

While the self-similar properties of the Jacobians of those conformal mappings appear 
clearly in Fig.^], one could wonder how the local symmetry breaking affects the continuous 
problem. We can see that if D~^{Q is univalued along the embedded tree, D:^^{C.) does vary, 
what makes the tree locally nonuniform and leads to a multifractal behavior. In other words, 
different paths of same length along the tree have the same weights in the symmetric case, but 
have different ones in the nonsymmetric case. The probability of a random path C of length L 
can be written in terms of a path integral with a Wiener measure 
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PC = V{s}exp 



D[s{t)] 




dt 



(47) 



where s{t) is a parametric representation of the path C. 

The first horocycles in Fig.^] can be parameterized as follows 

= i±(i-f sin0) 
= ^(l-cos^) 



(48) 



with 6 running in the interval [0,7r/2]. The condition ensuring the constant velocity s = 
along the horocycles gives with ( PD 

1 de 

— = const 

r] dt 



ds 
dt 



hence 



9{t) = arctan ( - 



(49) 



with proper choice of the time unit. This parameterization is used to check that the embedded 
tree is isometric. Indeed, the horocycles shown in Fig.^ correspond to a periodic sequence of 
steps like /3i/32/3i/32..., P1P3P1P3... or /32/33/32/93.... It is natural to assert that a step carries a 
Boltzmann weight characterized by the corresponding local values of D^. Therefore the period 
of the plot shown in Fig.^ is directly linked to the spacing of the tree embedded in the profile 
D-}. 



l/D 



ns 




. 65 



6 



FIG. 9. D^^ running along a horocycle at constant velocity. Periodicity shows that the tree is 
isometric. 
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Coming back to the probability of different paths covered at constant velocity, one can write 

'•ia dt 



log PC oc / 



D[s{t)] 



(50) 



The figure |T0| shows the value — logpc in symmetric and nonsymmetric cases for different paths 
starting at ti = 0^ and ending at t. In the symmetric case all plots are the same (solid line), 
whereas in the nonsymmetric case they are different: dashed and dot-dashed curves display 
the corresponding plots for the sequences /32/33/?2/?3--- and PiP^PiPs.... 

Following the outline of construction of the fractal dimensions Dg in Section ||, we can 
describe multifractality in the continuous case by 



1 1 
^""^ InAr(L) 



In 



T>{s} exp 



-q 1^"^ D~,\s{t))dt^ 



Pjs} exp 



L 



D^\s{t))dt 



(51) 



where A/'(Z) is the area of the surface covered by the trajectories of length L. This form 
is consistent with definitions |1| and |T6[ Indeed, if instead of the usual Wiener measure one 
chooses a discrete measure dxr, which is nonzero only for trajectories along the Cayley tree, 
we recover the following description. 

Define the distribution function 6(/5i, /?2, /?3, ^) = ®(^'^'^)' "which has sense of the 
weighted number of directed paths of k steps on the nonsymmetric 3-branching Cayley tree 
shown in Fig.^ The values of the effective Boltzmann weights |^ and || are defined in terms 
of the local heights of the surface along the corresponding branches of the embedded tree. 
We set 



Pi 
Ps 



— = exp 



exp 



t2 



dt 



Jti Dns[Srit)] 

*2 dt 



dt 



t, D^,[si{t)] Jt, D„s[s/(t)] 



t2 Dns[Srit)] 
*3 dt 



i 1.07; 
1.19 



(52) 



where ti,t2, ts are adjusted so that Sr(t) represents a step weighted with /^s for ti < t < t2 and a 
step weighted with /5i for t2 <t <t^ for right-hand-side horocycles while si{t) represents a step 
weighted with /?3 for ti < t < t2 and a step weighted with P2 for t2 < t < t-s for left-hand-side 
horocycles. 
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FIG. 10. Probability of different paths along horocycles. Dashed and dot-dashed curves correspond 
to right-hand- and left-hand-side horocycles on the nonsymmetric surface, while solid line shows both 
right-hand- and left-hand-side horocycles on the symmetric surface. 

The partition function B can be computed via straightforward generalization of 

Eq.@; it can be written in the form: 



(53) 



where Ai, A2 and A3 are the roots of the cubic equation 

A^' - A f 1 + ^ + 



PI , Pll32\ (Pi 132 , 



2 +^ =0 



V PI PI ) \ Pi ■ Pz, 

and v4o,-Bo and Cq are the solutions of the following system of linear equations 
Ao + i?o + Co = l + ^ + ^ 

A , J2 , J1P2 



Pz p3 Pi 



^oAi + -B0A2 + C0A3 



Pi . P2 . PI , P! , ^PlP2 , P!P2 , PlPl 

Ps Ps Pi Pi Pi PI PI 



Knowing the distribution function O (1^, A;) , Eq.(|l|) with the discrete measure dxr reads 
now (compare to (p^)-([T7|)) 



1 InG ( 
lim ^ 

g — 1 k-^oo 



'Pi 


<? 


'P2 




.ft. 


5 


.ft. 



,A:)-glne(|,|,fc 



ln(3 X 2'=-!) 

The plot of the function Dq{q) is shown in Fig.|TT| (the plot is drawn for k = 100 000). 



(54) 
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FIG. 11. Multifractality of trajectories on nonsymmetric tree with Boltzmann weights defined by 
the Jacobian of the conformal mapping 



V. DISCUSSION 



The results presented in Sections p|^|V| are summarized; they underhne several problems 
still unsolved related to our work, and raise the issue of their possible applications to real 
physical systems. 

1. The basic concepts of multifractality have been clearly formulated mainly for abstract 
systems in |1[]. In the present work, we have tried to remain as close as possible to these 
classical formulations, while adding to abstract models of Ref. [Q the new physical content of 
topological properties of random walks entangled with an array of obstacles. Our results point 
out two conditions which generate multifractality for any physical system: (i) an exponentially 
growing number of states, i.e. " hyperbolicity" of the phase space, and (ii) the breaking of a local 
symmetry of the phase space (while on large scales the phase space could remain isotropic). 

In Section |I| we have considered the topological properties of the discrete "random walk 
in a rectangular lattice of obstacles" model. Generalizing an approach developed earlier (see 
for example [0 and references therein) we have shown that the topological phase space of the 
model is a Cayley tree whose associated transition probabilities are nonsymmetric. Transition 
probabilities have been computed from the basic characteristics of a free random walk within 
the elementary cell of the lattice of obstacles. The family of generalized Hausdorff dimensions 
Dq{q) for the partition function f2(/3'^, k) (where k is the distance on the Cayley graph which 
parameterizes the topological state of the trajectory) exhibits nontrivial dependence on g, what 
means that different moments of the partition function f2(/3'', k) scale in different ways, e.g. 
that f^(/3'^, k) is multifractal. 

The main topologically-probabilistic issues concerning the distribution of random walks in 
a rectangular lattice of obstacles have been considered in Section |I|. In particular we have 
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computed the average "degree of entanglement" of a A^-step random walk and the probability 
for a A^-step random walk to be closed and unentangled. Results have been achieved through 
a renormalization group technique on a nonsymmetric Cayley tree. The renormalization pro- 
cedure has allowed us to overcome one major difficulty: in spite of a locally broken spherical 
symmetry, we have mapped our problem to a symmetric random walk on a tree of effective 
branching number z depending on the lattice parameters. To validate our procedure, we have 
compared the return probabilities obtained via our RG-approach with the exact result of P.Gerl 
and W.Woess [|15| and found a very good numerical agreement. 



The problem tackled in Section |I^ is closely related to the one discussed in Section |I[ We 
believe that the approach developed in Section |V| could be very important and informative as 
it explicitly shows that multifractality is not attached to particular properties of a statistical 
system (like random walks in our case) but deals directly with metric properties of the topo- 
logical phase space. As we have already pointed out, the required conformal transforms are 
known only for triangular lattices, what restricts our study. However we explicitly showed that 
the transform z^J^Q maps the multi-punctured complex plane z onto the so-called "topologi- 
cal phase space" , which is the complex plane C, free of topological obstacles (all obstacles are 
mapped onto the real axis). We have connected multifractality to the multi-valley structure 
of the properly normalized Jacobian Das{C,,i]) of the nonsymmetric conformal mapping Zns{C)- 
The conformal mapping obtained has deep relations with number theory, which we are going 
to discuss in a forthcoming publication. 

2. The "Random Walk in an Array of Obstacles" -model can be considered as a basis of a 
mean-field-like approach to the problem of entropy calculations in sets of strongly entangled 
fiuctuating polymer chains. Namely, we choose a test chain, specify its topological state and 
assume that the lattice of obstacles models the effect of entanglements with the surrounding 
chains (the "background"). Changing Cx and Cy one can mimic the affine deformation of the 
background. Investigating the free energy of the test chain entangled with the deformed media 



is an important step towards understanding high-elasticity of polymeric rubbers |Tl|. 

Neglecting the fiuctuations of the background as well as the topological constraints which 
the test chain produces by itself, leads to information losses about the correlations between the 
test chain and the background. Yet, even in this simplest case we obtain nontrivial statistical 
results concerning the test chain topologically interacting with the background. 

The first attempts to go beyond the mean-field approximation of RWAO-model and to 
develop a microscopic approach to statistics of mutually entangled chain-like objects have been 
undertaken recently in [Q. We believe that investigating multifractality of such systems is 
worth attention. 
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